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Abstract. Whereas the first part of this paper dealt with the relaxation in the /3-regime, this part inves- 
tigates the final relaxation (a-relaxation) of a simulated polymer melt consisting of short non-entangled 
chains in the supercooled state above the critical temperature T c of ideal mode-coupling theory (MCT). 
The temperature range covers the onset of a two-step relaxation behaviour down to a temperature merely 
q ■ 2% above T c . We monitor the incoherent intermediate scattering function as well as the coherent inter- 

' mediate scattering function of both a single chain and the melt over a wide range of wave numbers q. 

Upon approaching T c the coherent a-relaxation time of the melt increases strongly close to the maximum 
9max of the collective static structure factor S q and roughly follows the shape of S q for q > q max . For 
smaller g-values corresponding to the radius of gyration the relaxation time exhibits another maximum. 
The temperature dependence of the relaxation times is well described by a power law with a (/-dependent 
exponent in an intermediate temperature range. Deviations are found very close to and far above T c , the 
onset of which depends on q. The time-temperature superposition principle of MCT is clearly borne out 
in the whole range of reciprocal vectors. An analysis of the Q-decay by the Kohlrausch- Williams- Watts 
■ (KWW) function reveals that the collective KWW-stretching exponent and KWW-relaxation time show a 

modulation with S q . Futhermore, both incoherent and coherent KWW-times approach the large-q predic- 
tion of MCT already for q > g max . At small q, a q~ 3 -power law is found for the coherent chain KWW-times 
similar to that of recent experiments. 

PACS. 64.70.Pf Glass transitions - 61.25.Hq Macromolecular and polymer solutions; polymer melts; 
swelling - 61.20.Ja Computer simulation of liquid structure 
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© ■ 1 Introduction region where the dynamics becomes more and more dom- 

[ inated by activated processes over (free) energy barriers 

>> j. rfn , j c , , during supercooling (see pnotl for a detailed discussion of 

. -h -1 he preceding paper m reported results from a molecular- .. ° . , . » oi_ai 

k>( j • , c uj • jif this physical picture), 

dynamics simulation for a bead-spring model of a super- r J 1 ' 

5h 1 cooled polymer melt. The aim of this work was to extend An important property of T c is that it is an equilib- 

our previous analysis |,|,|,|,|,|| to coherent scattering rium quantity which can be calculated if accurate data for 

of the chains and the melt. Much of this analysis was the static structure factor are available. Such a quantita- 

guided by the mode-coupling approach to the structural tive approachwas performed for hard-sphere-like colloidal 
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glass transition |,0M|U). ' particles |2|,|lJ]^pft spheres |J or binary Lennard- Jones 
Mode-coupling theory (MCT) derives an equation of ( LJ ) mixtures @. These approaches reveal that MCT 
motion for density fluctuations (incoherent and coherent Provides a reasonable framework for quantitative predic- 
scattering functions), in which a non-linear coupling be- tl0ns of specific systems, but it also tends tooverestimate 
tween different wave vectors ( "modes" ) gives rise to a crit- the fr eezmg ability of a glass former (see |1§ for a com- 
ical temperature T c . This temperature is situated in the parative discussion). 

region of the supercooled liquid above the glass transition This tendency is exemplified by the so-called ideal- 
temperature T g . It may be interpreted as a crossover point ized MCT which deals with the cage effect only 



from the high-temperature region of structural relaxation • The idealized theory predicts the intermittence of 
(dominated by the "cage effect") to the low-temperature complete structural relaxation if T < T c . However, since 

T c > T g , the arrest of structural relaxation is in general 
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lations |8|,[18|]. Additional relaxation channels, which are 
ignored by the idealized MCT, must become important 
close to and particularly below T c . The extended MCT 
tries to include these channels approximately [^oHl9| ,^0|. 
These additional relaxation processes, called "activated 
hopping processes", replace the complete freezing of the 
idealized theory at T c by the aforementioned crossover to 
low-temperature activated dynamics. The status of this 
extension is, however, unclear at present On the 




one hand, evidence from simulations is scarce |B 2 
the other hand, evidence from experiments |2G 
challenged due to experimental problems [ p6[ . 

Nonetheless, a recent application of a schematic MCT 
model, including hopping processes, to light-scattering, di- 
electric relaxation and neutron scattering data [ pl[ yields 
a coherent description of the different spectra. This study 
also exemplifies the theoretical prediction [n]] that there is 
a temperature interval above T c , where the idealized the- 
ory represents a viable approach. The same conclusion was 
also drawn from various other experiment s [5 13^ 2^,|28l,p9| , 
I^HHD and simulations 




This has led to extension of the idealized theory to in- 
clude corrections to the asymptotic behavior close to T c 
J4l|,ff2| or to orientational degrees of freedom 
but also to further tests of the theory. 

In this spirit, we want to complement our previous 
simulation studies by an analysis of the a-relaxation of 
the coherent scattering functions of a supercooled poly- 
mer melt. This analysis is based on a detailed investiga- 
tion of the /3-relaxation described in the preceding paper 

The present paper is organized as follows: Section || 
introduces the model and the analysed quantities, whereas 
section || reviews the theoretical background for the sub- 
sequent analysis which is described in Section ^. The final 
section contains our conclusions. 



2 Model and Analysed Quantities 

The model we simulated by means of molecular-dynamics 
simulation is described in part I of this paper El (hereafter 
referred to as part I) and more extensively in |2j . Here, we 
only give a brief summary. 

We simulated a bead-spring model of linear polymers 
with N = 10 monomers per chain. The monomer-monomer 
interaction was modelled by a truncated and shifted Len- 
nard-Jones (LJ) potential. So, all quantities are measured 
in LJ-units in the following. In addition, bonds along the 
polymer backbone were introduced by a FENE-potential. 
The parameters of the potentials were adjusted such that 
the minima of the LJ-potential and of the bond poten- 
tial (= combined LJ- and FENE-potentials), r m j n and 
''bond, respectively, were geometrically incompatible. This 
means that the preferred nearest-neighbour distance of 
non-bonded monomers r m j n and the preferred bond length 
'"bond do not allow arrangement of the beads in a regular 
fee (or bec) lattice structure. Thus, cristallisation is pre- 
vented by this incompatibility and by the flexibility of the 
polymer backbone at all temperatures (the (square) end- 



to-end distance = 12.3 ±0.1 and the radius of gyration 
i?g = 2.09 ± 0.01 stay almost constant §). 

Simulations were done at constant pressure p = 1 
and at temperatures T = 0.46, 0.47, 0.48, 0.50, 0.52, 
0.55, 0.6, 0.65, 0.7, covering the range from the onset 
of glassy behaviour to about the critical temperature of 
mode-coupling theory (MCT), T c = 0.450 ± 0.005 (see §j 
and section 4.1). For all temperatures an amorphous static 
structure is observed, as discussed in part I. 

Whereas the first part of the paper investigated the 
behaviour of our system in the /3-relaxation regime, this 
part discusses the a-relaxation, i.e., the late-time struc- 
tural relaxation of the polymer melt. We computed the 
incoherent, the coherent chain and the coherent melt in- 
termediate scattering functions to probe the dynamics. 

We define the coherent (or collective) intermediate scat- 
tering function of the melt by 



<j> q (t) :-- 



s q (t) 

S q (0) 



(1) 



where S q (t) is given by 



M M 



S S) : = i(EE ex p{ i( i' Ht)-T S (p)]}) . (2) 

i=l j=l 

Monomer i has coordinate rj(t) at time t, M is the to- 
tal number of monomers in the melt, and (•) denotes the 
canonical ensemble average. S q (0) = S q is the melt's static 
structure factor. For defining the coherent chain scatter- 
ing function (i) we have to use N instead of M as upper 
summation index in equation (^) and sum over all pairs 
of monomers in a polymer. If we set i = j in equation (||), 
we obtain the single particle incoherent scattering function 
(j> s q (t). A superscript "s" for incoherent (self) and a super- 
script "p" for chain (polymer) quantities will be used in 
the following, while coherent quantities do not carry a su- 
perscript. When referring to all scattering functions, we 
shall frequently use the notation 4> q (t). The superscript 
"x" denotes the dependence on the correlators and is used 
analogously for other quantities. Further details about the 
computation of the scattering function may be found in 
part I. In our isotropic system all quantities depend on 
the modulus q — |q| of the reciprocal vector q only. 



3 Theoretical Background 

This section briefly summarizes predictions of the ideal- 
ized mode-coupling theory (MCT) for the a-process (see 
e.g. ||,[ll],[l8)), which we will apply in the analysis of the 
simulation data. The relevant temperature scale of the 
theory is defined by the separation parameter a, 



T,-T 



(3) 



Here, C is a system dependent constant. The separation 
parameter determines the temperature dependence of the 
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universal a-time scale by 



to 
\<t\i 



T>T C , 



(4) 



where to is the microscopic matching time of MCT. For 
our system, the critical temperature T c and the exponent 
7 of the a-process are given by T c = 0.450 ± 0.005 and 
7 = 2.09 ±0.07 §. 

Mode-coupling theory predicts that the correlator (t) 
decays in two steps in the vicinity of T c . In the first step 
<f>q(t) relaxes towards a plateau value, the non-ergodicity 
parameter f q c (0 < f q c < 1), whereas the second step 
represents the first part of the final structural relaxation 
to (a-process). 

An important prediction for the a-process is the time- 
temperature superposition principle (TTSP). This means 
that correlators (fr£(t), measured at different temperatures, 
fall onto a master curve in the limit a — > when time is 
rescaled by the a-relaxation time r q , defined by the decay 
of <fi q (t) to a certain value d, 



^(r 9 x ) = d . 



(5) 



The TTSP implies that the long-time relaxation of </)*(£) 



may be expressed as 



(6) 



Furthermore, the a-relaxation time r* is predicted to fac- 
torize into two parts: into a correlator (x) and wave vector 
(q) dependent part C x (g, d) and into the universal a-time 
f, 

to 



^(T) = C^(q,d)f(T) 



c*( g ,dy 



(7) 



Thus, all relaxation times r q in the a-regime exhibit the 
same temperature dependence. This property is in fact 
the reason why r* may be defined by 4>(t*) — d, with u d" 
being any arbitrary number sufficiently smaller than the 
non-ergodicity parameter f q c so that <fi{ T q) probes the a- 
process. In our analysis, we chose d — 0.1, which satisfies 
this condition well for all g-values under investigation. A 
different choice of d would only alter C*(q, d), but not the 
temperature dependence. 

The short-time expansion of the correlators in the a- 
regime is given by |4l| 



<f> x q (t) = f? - KB (t/f) 



(8) 



with the von Schweidler exponent b = 0.75 ±0.04 and B = 
0.476 ±0.060 j| for our system. The critical amplitude h q 
and f q c are investigated in detail in part I. 

In experiments and simulations one often finds that a 
Kohlrausch- Williams- Watts (KWW) function, defined by 



*f (*) = If 



exp 



(9) 



gives a very good description of the correlator in the a- 
regime. However, especially in the crossover region from 
the ft- to the a-regime small systematic deviations are 



commonly reported |ll| , |46| , |47 >|4g[|. Mode-coupling theory 
rationalizes these deviations by the fact that the short- 
time expansion of equation (0) does not coincide with 
equation (||) because b ^ ft q in general p8| . However, 
in the limit of large q it was proved that all correla- 
tors exhibit the following KWW behaviour, 

fim t%($ = /* c exp [-r*(t/f ) b ] ,r*^q. (10) 

So, / xK = / xc in this limit. On the other hand, the short- 
time limit of equation (@), which is valid for any q, suggests 
that the KWW-prefactor should always be equal for at 
least close) to / xc . Another consequence of Eq. (|lfj) is 
that the stretching exponent /3 xK satisfies 



lim I3* K = b 



(11) 



and that the g-dependence of the KWW-time-scales t xK 
is given by 



lim r xK = (r*y 1/b T otq-^f 



(12) 



The predictions presented in this section are valid in lead- 
ing order of er, corrections are of order |er|. This is the 
reason why the a-scaling should extend to higher temper- 
atures than the /3-scaling, corrections to which are of order 

§[![§■ 

On the other hand, deviations are in general expected 
at high temperatures, where a is large, and at tempera- 
tures close to T c , where relaxation processes, which are 
not taken into account in the idealized MCT, begin to 
dominate the dynamics [^9[. 



4 Results 



4.1 a-Relaxation Times 



The a-relaxation times r x can be read off from the simu- 
lated scattering functions by determining the time when 
^</( T <?) = 0-1 ( sce Eq. (]|)). This means that the a-decay 
times are obtained independently of any other quanti- 



•xK 



Figure [l] shows the collective 



ties, such as or ft 
a-relaxation time r q and the static structure factor S q 
as a function of q. The comparison of both quantities re- 
veals two conspicuous features: First, as temperature ap- 
proaches T c , r q develops modulations which closely resem- 
ble the oscillations of S q if q > q m ax- This is particularly 
prominent at the maximum of the structure factor and 
also visible around the first minimum of S q . It means that 
the local environment around a particle determines the 
relaxation more and more as temperature approaches T c . 
Such an adjustment of the (/-dependence of r q to that of 
S q for T — > T c is predicted by mode-coupling theory for 
hard s phe res f4§| and also observed in other simulations 
[ ^3|[50|,pl|l and experiments on non-polymeric glass form- 
ers ||52|,|53||. In this respect, the considered polymer melt 
behaves similarly to simple liquids. Moreover, r q exhibits 



4 



M. Aichele, J. Baschnagel: a-Dynamics of Glassy Polymer Melts 



2ji/R 



co 



; } 


1 1 1 

4 




m 




J *■ 



_ S q (T=0.46) 
++x /1000, T=0.46 
M T q /100, T=0.52 
♦^T q /10,T=0.7 



15 



Static structure factor S q and collective a-relaxation 



Fig. 1 

times r q (T) (scaled) for T = 0.46, T = 0.52, and T = 0.7 
versus q. Vertical long dashed lines are drawn at 2ir/R g = 4.35 
and at g max ~ 7.15, the maximum of S q . S q does not change 
much in this temperature interval (see part I). So, it is shown 
for T = 0.46 only. 



an additional characteristic feature. For all temperatures 
there is a broad peak around q « 4.5, which can be in- 
terpreted as corresponding to the radius of gyration R g 
(reciprocal length 2ir/R g — 4.35). A similar peak at a 
wave vector corresponding to the size of a molecule may 
also be observed in a simulation of a diatomic liquid |33| . 

Although the abovementioned features are clearly vis- 
ible in our simulation, we are not aware of similar findings 
in experiments on glass forming polymer melts . W e want 
to postpone a discussion of this point to Section 4.3, where 
we show results on KWW-fits. The KWW-function is also 



used for the description of the cv-decay in experiments |p4j 
|55j so that a qualitative comparison is better possible than 

for Tq. 

Figure || tests the prediction of equation (fjj) for the 
a-relaxation times of the coherent scattering functions. It 
shows tP and r q versus T—T c in a double logarithmic plot, 
where T c is kept fixed at T c = 0.450. In this representa- 
tion, MCT predicts the simulation data to fall on straight 
lines with slope 7 = 2.09 for all relaxation times, irrespec- 
tive of the g-values and the correlator considered. This 
expectation is not fully borne out. The data are linearized 
in an intermediate temperature interval only. Deviations 
appear for all q- values if T — T c < 0.02, and for large wave 
vectors if T — T c > 0.1. These deviations are not unusual. 
They are also observed in other simulation studies |33| , p0| 
and can be rationalized as follows: Equation (R) is only 
valid asymptotically close to the critical point so that de- 
viations must occur as T — T c grows. The onset of these 
deviations depends on the correlator and on q. Large q 
probe small distances in real space where microscopic de- 
tails should be felt more strongly. So, it is possible that the 
temperature region of applicability is left earlier at large 
q. On the other hand, very close to T c the idealized MCT 
is no longer expected to be generally applicable because it 
neglects certain relaxation channels (called "hopping pro- 
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Fig. 2. a-relaxation times of the chain, r q (top), and the 
melt, T q (bottom), versus T — T c in a double logarithmic plot 
for selected q-values as indicated. Symbol sizes are one stan- 
dard deviation. In this plot, T c = 0.450 is kept fixed. Accord- 
ing to equation (fjj), the plot should yield straight lines with a 
slope 7 = 2.09. Furthermore, MCT expects the power law to 
extend to larger temperatures than those, where the /3-process 
can be observed UM. In accord with the latter expectation, 
the simulation data exhibit a power-law behaviour for T < 0.7 
(except for the largest q- values), whereas a quantitative anal- 
ysis of the /3-process was only possible for T < 0.52 (see part I 
and Q). However, contrary to the MCT prediction, there is 
a non-negligible variation of 7 with q. The solid lines show 
least-square fits, including the maximum number of temper- 
atures, to determine the exponent 7. The dashed line shows 
(T - T c )" 2 ' 09 for comparison. 



cesses @M[56[|7|]), which are present in (almost) all 
real systems (some colloidal suspensions being an excep- 
tion; see [^3[|2^] or for a recent review). These pro- 
cesses allow the correlators to decay even below T c , where 
the ideal MCT would have predicted a cease of complete 
structural relaxation 

The region where the relaxation times exhibit a linear 
behaviour in Figure || were used to calculate the exponent 
7 with least square fits. The same analysis had already 
been done for the incoherent scattering function in ||. 
Figure || displays the combined results as a function of 
the wave numbers q. Qualitatively, the same behaviour is 
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observed for all correlators. There is a systematic drift of 
7 with q, which is not accounted for by (the leading-order 
redictions of) ideal MCT. The outcome of the /^-analysis 
7 = 2.09, only represents a good approximation for 
3 % Q ;$ 8. For large q it underestimates the adjusted 
7- values, whereas it overestimates them when q — ► 0. In 
this low-q limit the simulation results tend towards the 
exponent jo determined from the diffusion coefficient of 
a chain (D ~ (T — T C )^ D ) , i.e., 7D = 1.84 ±0.02. 

To cross-check the results of Figure || we also plotted 
(t^)" 1 / 7 versus T with fixed 7 = 2.09. From equation (^) 
it is expected to find straight lines for each q, which in- 
tersect the temperature axis at T c . As in Figure || visible 
deviations from a power law occur for all q at T — 0.46 
and at large q for the highest temperatures. Applying least 
square fits to the data points for 3 < q < 16 in the temper- 
ature interval 0.47 < T < 0.52, for which the /3-analysis 
was possible (see (l)), we found a systematic shift of T c to 
smaller values with decreasing q. Qualitatively, this find- 
ing resembles that of 7. However, the g-dependence of the 
critical temperature is weaker than that of 7 so that the 
results are compatible with T c = 0.450 ± 0.005 within the 
error bars if q > 3. This agrees with the analysis of the 
incoherent scattering function |||,0|. 



4.2 Time- Temperature Superposition Principle 

Figure ^ tests the validity of the time-temperature su- 
perposition principle (TTSP) expressed by equation @. 
It shows the coherent scattering functions versus time 
rescaled by the respective a-relaxation times, i.e. <p q (t/T q ), 
for q = 3.0, 6.9, and 14.0 at all investigated temperatures. 
The relaxation times are defined by <fi q {T q ) = 0.1. For this 
choice all correlators, even at high q, are sufficiently in 
the a-regime so that the /3-process does not perturb the 
behaviour. In the late ct-regime the curves collapse onto a 
master curve, with visible deviations for q = 14.0 at the 
highest temperatures. Figure || already showed that the 
power law for the relaxation times is violated for q = 16 
at T > 0.65. Thus, it is not surprising to find deviations 
from the TTSP master curves for similar q and T in Fig- 
ure ||, too. 

With decreasing temperature the TTSP property starts 
at earlier rescaled times and extends well above fjj 
suit often obtained from computer simulations | 

This trend is vio 



atec 



for 



and numerical calculations 
the lowest simulated temperature T = 0.46: Especially in 
the early a-region the curves for T = 0.46 deviate from 
the master curve. We believe that this observation is not a 
non-equilibrium effect as we compared <j> q (t) from 10 sim- 
ulation runs shifted in time by about 10 x T 9max and found 
no significant differences. An explanation could rather be 
that equation (||) is no longer valid because T = 0.46 is so 
close to T c that relaxation processes, which are not taken 
into account by idealized MCT, begin to dominate the 
dynamics. For the /3-relaxation the influence of these ad- 
ditional processes (called "hopping processes" by MCT) 
have been worked out in the extended mode-coupling the- 
ory [091. This analysis shows that these relaxation pro- 
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Fig. 3. Q-relaxation exponents 7 from the incoherent scat- 
tering function (top, from and the coherent counterparts 
versus q extracted by least-square fits (see Figure g). 7 = 
2.09(±0.07), the result from the /3-analysis ^j, is shown as 
a thick solid line. To indicate the q — > limit the exponent, 
70 = 1.84(±0.02), from an analysis of the diffusion coefficient 
of a chain is drawn as a thick dashed horizontal line. 



cesses change the shape of 4%(t) in the late-/3/early-a 
regime and remove the power-law divergence of the a-time 
scale f (see Eq. Therefore, it is possible that devia- 
tions from the idealized TTSP occur if T is very close to 
T c . Such deviations have already been observed for <jf q {t) 
for the present model B. 

For q < 4 an oscillation at the beginning of the 0- 
regime can be observed at all temperatures, which be- 
comes more pronounced at lower T. With decreasing tem- 
perature the amplitude of the oscillation increases when q 
becomes smaller, as exemplified in Figure for T = 0.47. 
Sound waves travelling through the simulation box and 
reentering again at the opposite side could be the reason 
for this effect j58|. The sound velocity was estimated to 
be c w 7 (§. Using 2n/cq = t° sc we obtain times t° sc 
comparable with the times where the oscillations occur in 
Figure ||. Note that neither cjf q (t) ||) nor </)? (t) show such 
oscillations. This might be a hint that a non-linear cou- 
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Fig. 4. Demonstration of the time-temperature superposi- 
tion principle (TTSP) for the coherent intermediate scatter- 
ing function 4> q (t) for all investigated temperatures, T = 0.46, 
0.47, 0.48, 0.50, 0.52, 0.55, 0.6, 0.65, 0.7 (from left to right in 
the plots) for q — 3.0 (top), q = 6.9 (middle), and q — 14.0 
(bottom). At T = 0.46 (dot-dashed line) the TTSP does not 
work as well as for the other temperatures. The a-relaxation 
times were taken to be the decay time to 0.1 (dotted horizontal 
line). For each q the non-ergodicity parameters are shown 
as dashed lines. 
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Fig. 5. Coherent intermediate scattering function (f> q (i) at 
T — 0.47 for the lowest investigated wave-numbers q up to 
Qmax = 6.9. At small q an oscillation at 0.15 < t < 1.0 develops, 
which is seen neither for <j> q (t) 0] nor for <fi q (i). 



pling between different q- values of the coherent scattering 
function is important for their presence. 

Another feature of the data are weak bumps at the 
beginning of the /3-plateau, which occur at t « 3 for all 
correlators and all q if T < 0.48. These arc discussed in 
[Q for <p s q {t). They are hard to see in Figures |^ and || for 
T — 0.47 because the statistics is much worse for 4> q (t) 
than for 4*„(t). So, it is difficult to distuinguish them from 
statistical fluctuations. In part I a weak bump is also vis- 
ible in Figures 4 and 12 at t w 2. As the aforementioned 
oscillations these bumps could be related to sound waves 
[HsJ . Using L = 10.5 for the linear dimension of the sim- 
ulation box and c ~ 7, one obtains t — L/c « 1.5, which 
is close to the position of the bumps. If they were really 
related to sound waves, their position should shift when 
changing the size of the simulation box. We have not at- 
tempted such a finite size study yet, although simulation 
results on two-dimensional polymer melts ]59 ], on binary 
mixtures of LJ-particles |6(| or of hard spheres |6l[] and 
on silica glasses ^||,^3) suggest that there might be pro- 
nounced effects at low temperature. 

Besides sound waves the oscillations could still have 
other causes, for instance, non-equilibrium effects fMj . 
However, this cause seems unlikely, since the oscillations 
are equally present at all investigated temperatures and 
great care was taken to achieve equilibrium at the lowest 
temperatures. Furthermore, the Nose-Hoover thermostat 
used to maintain constant temperature in the simulation 
can also be ruled out as a crucial factor. In Ref. (|] its 
influence was thoroughly tested against simulations in the 
microcanonical ensemble. No significant difference of the 
dynamics was found. 

Of course, it is important to assess the relevance of the 
oscillations for the results in the a-regime. The analysis of 
the /3-process in part I showed that the oscillations perturb 
the /3-master curve only at the beginning of the /3-regime 
for q < 3 (see Figure 5 of Q). For longer times no de- 
viations can be observed in comparison to the behaviour 
of the correlators at larger q- values, where the oscillations 
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are absent. If the major part of the /3-process is already 
unaffected, it is unlikely that the oscillations could have a 
pronounced influence on the a-relaxation discussed here. 



4.3 Kohlrausch-Williams-Watts Fits 

When performing KWW-fits one has to keep in mind that 
the results will vary quite significantly with the chosen fit 
interval (see below and also |^,^5|). Since the /3-process 
overlaps with the a-process, one would have to use a com- 
bination of equations (|^) and (|^) in order to incorporate 
the influence of the /3-dynamics. To alleviate this prob- 
lem sometimes only the late a-process is fitted ||§[5(J to 
minimize interference with the /3-process (see also W for 
a more detailed discussion of this and other approaches). 
Apart from these fundamental considerations it is clear 
that a highly non-linear fit will sensitively depend on its 
input data. 

We fitted KWW-functions using two different proce- 
dures: One method was a fit with three adjustable param- 
eters (fg K , r* K , /3* K ). This produces a description which 
is not biased by any theoretically motivated restrictions 
on (some of) the fit parameters. In the other method f* K 
was taken to be the non-ergodicity parameter (NEP) f q c . 
This choice is suggested by the MCT prediction (||). It 
also fixes the KWW-times via 4^{ T q K ) — fq C / e so that 
the stretching exponent /3* K is the only free parameter. 
In the following the results of both approaches, which we 
call "free fit" and "fixed-NEP fit" , are compared. 

The constraint /* K = /™ is not a very strong restric- 
tion: the fitted values for /* K are within 10% of the non- 
ergodicity parameter /* c . The deviations of f q K relative 
to f*° depend on both q and the size of the fit interval 
(= ^cut, see below) . If q > q max , f q K is almost indepen- 
dent of the chosen fit interval and systematically larger 
than /* c . For q < g max the influence of the fit interval be- 
comes bigger. Furthermore, f* K < f q c becomes possible. 
Despite these (quantitative) differences both f q K and /* c 
exhibit the same qualitative g-dependence. Consequently, 
the KWW-times from free and fixed-NEP fits also devi- 
ate only quantitatively (by at most 15% for q > 4), but 
not qualitatively. Both fit procedures yield the same q- 
dependence for the relaxation time. We will return to this 
point below. 

For the fits all data points with (f>g(t) < x cut fg C were 
used. The fit interval was systematically varied by increas- 
ing the cutoff parameter a; cut from 0.3 to 0.9 in steps of 
0.1 to investigate the influence of the chosen interval on 
the resulting KWW-parameters. Removing data below a 
certain small positive value leads to systematically higher 
stretching exponents, especially at large q. Therefore, we 
did not introduce a lower bound, although noise around 
zero might influence the results. However, such a lower 
bound - data points below 0.02 were discarded - was used 
in the study of (ff q (t) in 0, resulting in greater values for 
j3g than in our analysis. Nevertheless, the difference is less 
than 5% at the largest q. So, the statements made in [Q re- 
main valid. Because of the TTSP the shape of a-relaxation 
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Fig. 6. KWW-exponent f3 q of the coherent scattering func- 
tion 4> q (t) at T — 0.47 for different cutoff values x cut fq. For a 
given 2: cu t all data points with <j> q < x cu tfq were used for the 
fits. A thick dashed line is drawn at b = 0.75, the prediction of 
MCT for q -» oo (Eq. 0)). In the upper panel all three KWW 
parameters were fitted. In the lower panel only f3 q was ad- 
justed, fg — fq was fixed and thus also the time constant 
by <t>q(T q ) = fq/ e - One sees that 0^ is qualitatively in phase 
with the structure factor S q , especially when all three free pa- 
rameters are allowed to vary freely. The line types for both 
panels are identical. In the upper one the curve for x cut = 0.3 
is missing because the fit algorithm did not converge. 



is identical for all temperatures 0.47 < T < 0.52 (cf. Fig- 
ure @). Thus, we confined the KWW-analysis to T = 0.47. 

A comparison of the results for the stretching exponent 
f3 q obtained from fits with three free parameters (top) 
and with fixed = /| (bottom) is shown in Figure |^. 
Qualitatively, both plots are comparable. They show that 
0^ is roughly in phase with the static structure factor, a 
result also found in calculations for hard sphere systems 
p8[ , simulations of water p lj o r neutron scattering ex- 
periments of orthoterphenyl |29|. On the other hand, the 
data indicate a further feature which is not present in S q . 
There is a weak shoulder at q w 4.5 w 2ir/R g . This seems 
to correspond to the large r q encountered in Figure ^. 

Let us now compare both fit methods in more detail. 
The resulting values of 0^ differ especially at small and 
large q. At small q, this is not surprising, since the differ- 
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ences between fitted /* and f q c are most pronounced 
in this g-regime. In the case of the fixed-NEP fit, 0^ 
increases monotonously with the cutoff parameter when 
Xcut > 0.6, whereas it does not significantly depend on it 
for a; cu t < 0.5. This is an indication that the presence of 
the /3-process starts to be felt at a; cu t > 0.6. As the other 
two parameters are fixed, a better adaption to the /3-decay 
can only be accomplished by an increasing stretching ex- 
ponent (note that b — 0.75 > 0£ for these small q- values). 
Contrary to that, the fits with three free parameters do not 
show a systematic dependence on x cu %. The curves scat- 
ter randomly by roughly 15%. The size of this variation is 
comparable to that of the fixed-NEP fit. Nevertheless, the 
sensitivity to x cu t and the different results of both meth- 
ods illustrate that the choice of the fit conditions can lead 
to (significantly) different results, although the statistical 
accuracy of the input data is good (see Fig. ^) . 

As q approaches g max , the /S^-values of both methods 
converge and essentially agree with each other at g max , 
where stretching is least pronounced. However, if q in- 
creases beyond g max , the results of both methods differ 
again substantially. The fixed-NEP fit exhibits a strong 
dependence on x cu t, which is quite generally much larger 
than that of the free fit, and gives results close to the 
large-g limit 0^ = b (see Eq. (0)), whereas the /3^-values 
of the free fit are distinctly smaller and do not seem to 
tend to b. We only see an increase of 0^ at the largest q. 
This might suggest either that the previously determined 
exponent parameter A is too small or that the predicted 
large-g limit is approached from below - if at all - in our 
system. The latter behavior would be in contrast to cal- 
culations for hard spheres (4f| and simulations of water 
|^l|,^6). A test of this conjecture would require an exten- 
sion of the analysis to greater q. However, this extension 
is hampered by the very small plateau values at large g, 
which renders it difficult to discriminate the a-decay from 
the statistical noise around zero. 

In calculations for ideal hard spheres the limit of 0^ = 
b is reached at q rj 6 x q max For our system this 
would correspond to q w 42. This might be the reason 
why /3q does not approach the large-g limit b in the range 
of wave vectors investigated for most fit intervals in the 
upper panel of Figure ||. Similarly, a KWW-analysis of the 
incoherent scattering of a binary Lennard-Jones mixture 
could not confirm the large-g behaviour of (3 q K by means 
of a free fit procedure ||(|. However, the lower panel of 
Figure || suggests that the choice 0f = b is not unreason- 
able if q > g max . This conclusion was also drawn from a 
fixed-NEP analysis of incoherent scattering in j| . There- 
fore, we performed the comparison of both methods again 
for <j)g(t) and also for coherent chain scattering. The re- 
sults for (3g K and (3^ K are the same as in the case of 0^. 
The behavior observed in Figure |^ is thus representative 
of all scattering functions of our model. 

Therefore the question arises of whether a KWW- func- 
tion with = f q and 0^ = b can provide a comparable 
description of the simulation data as a free fit at large g. 
Figure ^ shows a comparison of both approaches for 4> q (t) 
at various q > g max . At first sight, both methods yield 




t 



Fig. 7. Comparison between cf> q (t) from simulation data (sym- 
bols), KWW- fits with three free parameters to data points 
with 4> q (t) < fg/2 (dashed lines) and KWW-functions calcu- 
lated from the large-g prediction of MCT: lim q _>oo 0^ = b and 
4>(t^) = fq/e (solid lines) at some q > g max as indicated. The 
temperature is T = 0.47. t a is the /3-relaxation time [Q. The 
figure shows that the large-g approximation ( |lC)|) works quite 
well, but from the inset it can be seen that it does not describe 
the late a-process as well as the free fit. 

very satisfactory approximations for the a-decay, but \ 2 
is typically a factor of 10 smaller for the free fit. The origin 
of this difference is illustrated in the inset of the figure: 
Using ff = Pf = b, and from «^(r ? K ) = f c q /e, 
the description of the data works quite well for q > g max , 
except for the very late a-process which are better de- 
scribed by the free fit. The discrepancy does not diminish 
if /3* K is adjusted instead of posing 0f = b while still keep- 
ing = /£. The same results are also obtained for the 
incoherent scattering and the coherent chain scattering 
functions. So, we may conclude that the free fit procedure 
yields a better description of the simulation data than the 
fixed-NEP fit and that we cannot unambiguously confirm 
the asymptotic limit (3 q K = b in the g-range accessible to 
our simulation. 

However, the situation is different for the KWW-times. 
These times can be defined by (f>q( T q K ) — fq C / e ( or by 
0q( T q K ) — fq K / e )y independently of 0f. Since the q- 
dependence of f* c closely matches that of f* K and both 
KWW-prefactors are numerically not very different, the 
KWW-times from the free fits show the same behaviour 
as those of the fixed-NEP fits. So, only the KWW-times, 
defined by 0^(r^ K ) = fq C /e, are shown for T = 0.47 in 
Figure || and are compared with the collective a-relaxation 
times T q . If q < 12.8, r q and exhibit a qualitatively 
similar behaviour: There is a broad maximum around q m 
4.5 and a sharp maximum at g max followed by oscilla- 
tions which are in phase with each other and mimic those 
of S(q). However, the absolute values of both relaxation 
times differ. If g < 12.8, r q > t?, whereas r q < for 
larger q. At these large g-values the behaviour of both re- 
laxation times is also qualitatively different. r q decreases 
sharply, whereas all t* k curves fall onto a master curve if 
g > 11. This master curve is given by a power law which 
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Fig. 8. KWW-times r q K (symbols) defined by 
fq C /e from simulation data versus q in a double logarithmic 
plot. The q-values shown are: q — 2, 3, 4, 5, 6, 6.9, 7.15, 8, 9.5, 
11, 12.8, 14, 16, 19. r' K and r pK are accurate to within 5%, 



q < 12.8. On the other hand, if q > 12.8, the q-dependence 
of the right-hand side does no longer strongly oscillate, but 
is rather close to a power law which modifies that of t* k . 
This comparison shows that the asymptotic behavior of 
equation ( |l4| ) can only be observed if the q-dependence of 
the amplitude of the a-process is taken into account, as it 
was also done in the theoretical derivation of equation ( |l4| ) 

s 

At this point, we would like to note that, to the best 
of our knowledge, neither a modulation of with S q nor 
a prepeak at wave vectors corresponding to the radius of 
gyration has been reported in neutron scattering exper- 
iments. For instance, experiments on polybutadiene [ |54| 
or polyisobutylene [|55| rather find evidence for a power 
law q~ n (n = 3.1 |55| , n — 3.6 [|^|) in the range of the 
first and second maxima of S q . The power law is superim- 
posed with oscillations that seem to in phase with S q for 
polyisobutylene 55 , 67j , but in anti-phase for polybutadi- 
ene |Q. A possible origin of this difference between ex- 
periment and simulation could be the temperature range 



r 9 K to within 10% (two a intervals). The thick dashed line is studied. In the simulation, the modulation of T q (or r ( 
the large-g prediction of MCT, t* k cx q- 1/b (Eq. @), which 
is seen to be nicely fulfilled for all correlators. S q (scaled and 
shifted) is indicated with a dotted line. The a-decay times r q 
(defined by 4> q {j q ) = 0.1) are also shown as a dash-dotted line. 
Except at large q (q > 12.8) they exhibit the same qualitative 
behaviour as the KWW-times r q . Power laws q~ 2 and q~ 3 are 
shown for comparison. 



is very close to q 4 / 3 , the asymptotic result of MCT (see 
Eq. @). 

This disparity between r* and r^ K can be rationalized 
by the following argument. According to equation (^) the 
times 

r x = C x (g,d = 0.1)f; rf = C x (q, /»/e)f , (13) 

are related to one another. Assuming a KWW-behaviour 
with / xK = / xc for the correlators in the a-regime, we get 



r* _ C*(q,d) 
r xK C x (g,/-/e) 



IT > d 



(14) 

Above relation was tested for the coherent (chain and 
melt) a-decay and KWW-times. In the range 3 < q < 14 
and for all x cu t the deviation is smaller than 10%. This is 
a strong consistency test because the a-times t x were de- 
termined completely independently of the KWW-times, in 
which the non-ergodicity parameters from the /3-analysis 
also enter. If q < 12.8, the choice d = 0.1 implies that 
f q c /d > e so that r q > t^, and vice versa if q > 12.8. Fur- 
thermore, the previous discussion revealed that f q and p q 
are in phase with S q . By evaluating the right-hand side of 
Eq.(|l4) one sees that the ratio T g /r^ is modulated with S q 
as well. This modulation should be more pronounced for 
T q than for the KWW-times because the in-phase vari- 
ation of with Sn should amplify that of the right-hand 
side of equation (JlJ) . Figure || supports this expectation if 



with the coherent static structure factor only becomes vis- 
ible very close to T c . T = 0.47 corresponds to a reduced 
temperature of (T — T c )/T c w 0.044, whereas the measure- 
ments of S q of polybutadiene Q were done at 260 K, so 
at a reduced temperature of w 0.2 (using T C PB w 216 K 
[[38||6^]). In turn, this reduced temperature would corre- 
spond to T = 0.54 in our model, a temperature at which 
a modulation with S q is almost not present (see Figure Q) . 
The same conclusion should also hold for polyisobutylene 
because the temperatures investigated, T — 390, 470 K, 
are far above T g = 201 K J46fl . Furthermore, the prepeak 
on the scale of the radius of gyration has perhaps not been 
found in the experiments, since it would lie outside of the 
q-range studied in the coherent scattering of the melt. An 
estimate of the corresponding q- value with the experimen- 
tal value for R g (~ 19 A) for polyisobutylene suggests this 
conjecture. 

At small q, t| k increases less steeply than t^ k . This 
might be rationalized by taking into account that the 
chain has a structure, which decays in time, whereas no 
such additional relaxation happens for the single parti- 
cle correlator. At the smallest g-values (q < 5) rP K fol- 
lows approximately a q~ 3 power law behaviour, which 
is comparable to experimental results Free diffusion 
would require a q~ 2 behaviour, expected for q <§; l/i? g ~ 
0.7 [|[0|], which is not accessible in our investigation. At 
the lowest q-values, for which KWW-fitting was possi- 
ble, T q K increases more steeply than q~ 2 in agreement 
with previous findings 0]. Note that the investigated re- 
laxation times in 0] were defined by a decay to d = 
f q c /2 instead of d = f q c /e used here. According to equa- 
tion (fbil), this implies that the two times differ by a factor 
of C s {qJ q x /2)/C s {qJ s q c /e) « (ln2) 1 // 3 f. This function 
varies slowly with q so that the g-dependence of both re- 
laxation times should only be slightly different. 

Finally at q ss 6, r g K and t£ k cross over to the large-g 



prediction r x 



xK 



oc q 



-l/b 



(Eq. @) of MCT. Due to the 
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modulation of with S q the collective relaxation time 
of the melt enters the asymptotic regime later than r^ K 
or tP k (for q > 11). It is interesting to note that the 
large-q behaviour of MCT is seen rather clearly for all 
r* K , whereas the large-t? behaviour of /3* K could not be 
confirmed unambiguously - if at all (cf . Figure ^) . 



5 Conclusions 

This paper summarizes simulation results concerning the 
a-relaxation of a dense, non-entangled polymer melt. The 
melt was simulated by means of molecular-dynamics simu- 
lations, using a Lennard- Jones potential for modelling the 
monomer-monomer interactions and a FENE-potential for 
modelling the bonds. The simulations were done at tem- 
peratures 0.46 < T < 0.7 (measured in Lennard-Jones 
units), covering the temperature range from the onset of 
a two-step relaxation behaviour almost down to the criti- 
cal temperature T c ~ 0.45 of (ideal) mode-coupling theory 
(MCT). We computed the incoherent, coherent chain, and 
coherent melt intermediate scattering functions to probe 
both single-monomer and collective dynamics of the melt. 

The a-times t* and the KWW-times r* K are related 
to each other (see Eq. (|l4|)) and exhibit similar features 
when T approaches T c : a broad peak at wave vectors cor- 
responding to the chains' radius of gyration and a sharp 
peak at the maximum of the collective static structure fac- 
tor S q (Figures |] and ||). The peak at q w 2n/R g might 
perhaps be thought of as a consequence of the "packing" 
of the diffuse soft polymer coils in the melt, which should 
also contribute to the freezing of the system like that of 
the monomers does. This peak is absent in simple liquids. 
On the other hand, as in simple liquids, both time scales 
display a modulation that is in phase with S q for q > q max . 

The quantitative analysis of the temperature depen- 
dence of the coherent a-relaxation times shows that they 
follow a power law for T > 0.47 (Figure |). At T = 0.46, 
however, there are clear deviations. This implies that the 
temperature range, where idealized MCT is applicable, 
is left. Moreover, at large wave numbers the power-law 
behaviour is progressively violated at T > 0.65 when q 
increases. This temperature marks an upper bound of 
the applicability of the aymptotic predictions of idealized 
MCT. The results show that the bound is g-dependent. 
The exponent of the power law is only compatible with 
the MCT prediction from the /3-analysis for wave vectors 
close to q ma x, but in general grows with q for both inco- 
herent and coherent scattering functions (Figure ||) . This 
observation is (currently) not accounted for by MCT. 

On the other hand, the time-temperature superposi- 
tion principle, an important prediction of MCT, is well 
borne out by the simulation data if T > 0.47. For T = 
0.46, however, deviations from this scaling behaviour are 
encountered (Figure |4[). Within the framework of MCT 
such deviations are expected close to T c due to relaxation 
processes which are not taken into account in the ide- 
alized version of the theory. These processes change the 
shape of the correlators in the late-/3/early-a regime jl9| 



so that a superposition with curves measured at higher 
temperatures is no longer possible. However, MCT pre- 
dicts that the additional relaxation processes do not per- 
turb the scaling in the /3-regime. The factorization theo- 
rem should still hold [ pj| . Our simulation results provide 
evidence for this prediction (see Figure 5 of El). 

At the beginning of the /3-regime oscillations were ob- 
served for the collective scattering functions at all T which 
grow with decreasing q (Figures || and ||). Sound waves 
might be the origin of these oscillations, but a conclusive 
explanation is still to be found. It is important to note, 
however, that these oscillation do (presumably) not alter 
the findings for the a-relaxation presented here, since they 
already have no noticable influence on the major part of 
the (3- process (see part 10]). 

In a detailed analysis of the shape of the a-decay by 
the Kohlrausch- Williams- Watts function two procedures 
were compared: Firstly, we used three free fit parame- 
ters, and secondly, a fit with the KWW-prefactor set to 
the non-ergodicity parameter, thus fixing the KWW-time 
and only allowing the stretching exponent to vary. The in- 
dependence of the collective KWW-stretching exponents 
0^ shows a modulation with S q independent of the fit 
method (Figure ||). Furthermore, the asymptotic predic- 
tion of MCT for the q-dependence of the KWW-times, 
r xK oc is verified for all correlators for large q (Fig- 

ure ||). However, our data do not allow to clearly confirm 
the corresponding prediction for the stretching exponent 
in the same g-range. In case of the coherent chain scat- 
tering function, a q~ 3 dependence of the KWW-times at 
small q is obtained. This behavior is compatible with re- 
cent experiments [p5| . 

In summary, one can say that ideal MCT provides a 
good starting point for the description of the a-relaxation 
of our model if temperature is close, but not too close to 
T c . We could give estimates for the temperature range of 
validity of the theory and for the g-range where correc- 
tions to MCT predictions become important. There are, 
however, some polymer-specific effects, such as the strong 
increase of the a-relaxation time on the length scale of the 
radius of gyration, which cannot be contained in a theory 
developed for simple liquids. 

To what extent the results discussed in this and the 
preceding paper are also observable in simulations of real- 
istic polymer models is not clear at present. An important 
step in this direction was made by the development of a 
model potential for polybutadiene |71,7^|, which yields 
very g ood agreement with neutron scattering experiments 
[|73[|7j] . A comparison of the dynamics of this model in the 
supercooled state with our data would represent an impor- 
tant test of our findings. Results on a different model for 
polybutadiene f75| seem promising in this respect. There- 
fore, we hope that our simulations can contribute to devel- 
oping a detailed theory of the glassy behaviour of polymer 
melts. Attempts in this direction are underway [|76| . 
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